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Abstract 

Using Feynman path integral technique estimations of the ground state energy have 
been found for a conduction electron interacting with order parameter fluctuations 
near quantum critical points. In some cases only singular perturbation theory in the 
coupling constant emerges for the electron ground state energy. It is shown that an 
autolocalized state (quantum fluctuon) can be formed and its characteristics have 
been calculated depending on critical exponents for both weak and strong coupling 
regimes. The concept of fluctuon is considered also for the classical critical point 
(at finite temperatures) and the difference between quantum and classical cases 
has been investigated. It is shown that, whereas the quantum fluctuon energy is 
connected with a true boundary of the energy spectrum, for classical fluctuon it is 
just a saddle-point solution for the chemical potential in the exponential density of 
states fluctuation tail. 
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1 Introduction 



The physics of quantum critical point (QCP) [1,2,3,4,5] is now a subject of 
growing interest. There is a solid experimental evidence of relevance of the 
QCP and related phenomena for ferroelectrics [6], high-temperature super- 
conductors [7,8], Bose-Einstein condensed atoms in traps [9], itinerant elec- 
tron magnets [10,11,12,13], heavy fermion compounds [14,15] and many other 
systems. Similar to classical critical points or second-order phase transitions, 
a scaling concept is of crucial importance near the QCP and universal critical 
exponents can be introduced, which determine all anomalous properties of the 
systems near QCP [2]. The universality means that the basic physics depends 
not on the details of a microscopic Hamiltonian but rather on space dimen- 
sionality, dispersion law of low-frequency and long-wavelength fluctuations of 
an order parameter and symmetry properties of their effective action. In con- 
trast with classical phase transitions at finite temperatures thermodynamics 
of the QCP is essentially dependent on the dynamical critical exponents [1]. 

There is an interesting issue how these critical fluctuations can effect on the 
state of an excess charge carrier which appears as a result of doping, injection, 
photoexcitation, etc. One can consider for example the electron motion in a 
crystal near the ferroelectric quantum phase transition in virtual ferroelectrics 
such as SrTiOs or KTaOs under doping or pressure [6,16], or near quantum 
magnetic phase transition due to competing exchange interactions [2]. To our 
knowledge this problem has not been considered yet. One may speculate that 
a specific nature of the order parameter is not very essential for this prob- 
lem; due to softness and long-range character of the critical fluctuations the 
effects of their interaction with the conduction electrons may be very strong. 
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In particular, we will see that a self-trapping (autolocalization) of the carrier 
proves possible, similar to a polaron formation in ionic crystals [17,18] or spin 
polarons ("ferrons") in magnetic semiconductors [19]. A general concept of 
the self-trapped electronic state due to interaction with order parameter fluc- 
tuations ("fluctuon") has been proposed many years ago by Krivoglaz [20]. 
It appeared, however, that his phenomenological approach is not applicable 
near the critical point where the fluctuon radius is smaller than the correlation 
length [21]. We have considered this case [21,22,23] using Feynman path inte- 
gral variational approach developed him for the polaron problem [24,25]. Here 
we apply similar technique to consider the quantum case. It will be shown 
that the classical and quantum fluctuons are drastically different: if the latter 
can be considered as a specific quasiparticle the former one represents some 
quasilocalized state in the density of states tail. Apart from possible appli- 
cations to condensed matter physics the problem under consideration gives 
a nontrivial example of the interaction of a fermion with a bosonic quantum 
field with anomalous scaling properties. 

Whereas only the case of dispersionless Einstein phonon has been considered 
originally by Feynman, later this method has been used also to describe the 
interaction of electron with acoustic phonons [26,27]. We consider here a gen- 
eral case of fluctuations with arbitrary dynamics which can be, in particular, 
of dissipative type. The answers will be written in terms of some frequency 
momenta of the fluctuations. One can assume that the type of the fluctuation 
dynamics, being relevant, e.g., for transport phenomena is not essential for 
static characteristics such as autolocalization radius and energy; anyway, the 
method used by us gives a rigorous upper limit for the ground-state energy. 
Another difference (which is more important) is that the phonon field is Gaus- 
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sian whereas the Gaussian approximation for the fluctuations which we will 
use can be justified only for not too large coupling constants. It leads to some 
restrictions which will be derived separately for all cases under consideration. 

The interaction of electrons with quantum critical fluctuations are intensively 
studied, especially in connection with high-temperature superconductivity and 
heavy-fermion systems (for review, see Refs. [7,8] and [14,15], respectively). 
Usually it is assumed that the coupling constant is small in comparison with 
the Fermi energy. Here we consider the case of single carrier where the charac- 
ter of electron states is essentially different; one can say that this difference is 
similar to the difference between localized and extended states for the disor- 
dered systems. As a next step, it would be interesting to consider degenerate 
gas of fluctuons where the Fermi energy is finite but small in comparison with 
the autolocalization energy which might be a subject of future investigations. 

The paper is organized as follows. In Section 2 we overview general formalism 
for solving the problem posed. In the present paper we consider the case of not 
too large coupling constant, where the problem can be considered in Gaus- 
sian approximation for the interaction with the fluctuations; explicit criteria 
are presented below. The quantum case (zero temperature) is considered in 
Section 3. Using the scaling properties of the fluctuation spectral density (Sub- 
section 3.1) we construct regular perturbative expansion of the energy in the 
coupling constant (weak-coupling regime, Subsection 3.2) as well as singular 
perturbative expansion in strong coupling regime (Subsection 3.3). The very 
existence of the regular perturbative regime depends crucially on the value of 
dynamical critical exponent z and anomalous dimension d. The problem of 
fluctuon at classical critical point (finite temperature) is treated in Section 4. 
We solve the problem by both Feynman variational method (Subsection 4.1) 
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and using Green function technique with vertex corrections via Ward identity 
(Subsection 4.2). The similarity of the results as regards dependencies of the 
density of states on the energy and coupling constant justifies the variational 
approach. 



2 Formulation of the problem using Feynman path integral 

For simplicity, we will consider the case of a scalar order-parameter acting 
only on the orbital motion of the electron and not on its spin (for example 
it may be the QCP in ferroelectrics) . Then, in continuum approximation, the 
Hamiltonian of the system consisting of the electron and the order-parameter 
field can be written in a simple form 



H = H f M + He (r,tp) , H e (r,<p) = -\v 2 r - gy (r) (1) 

where we have chosen the units h = m = 1, to is the electron effective mass, 
r is the electron coordinate, <p (r) is the quantum order-parameter field with 
its own Hamiltonian 7if((p) and g is the coupling constant. The partition 
function of the whole system may be transformed to 



Z = Tre-WiM-PHeM = Z f (rr r T T exp - £ H e (r,y? (r,r)) dr 



) (2) 
/ 



where Zf = Tr^e P n f(^ is the partition function of the field, cp (r,r) 
e TH f^(p (r) e-^/W and 



(A(ip)) f = ^-1t v e-^A(ip) (3) 
z f 
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is the average over the field states. Using Feynman path-integral approach 
[25,28,29] and taking average over ip yields for the electron-only free energy 

T=-h\nZ-\nZ f ) = —\nf e s V[v(r)}, (4) 

P P ir(0)=r(/3) 

where So + Sint is the effective action, 



oo ,,0 ,p ,p 

Sint = -^—r ■■■ /C m (r(ri),ri;...;r(r m ),r m )dri...dr m (5) 
±r Q ml Jo Jo 



dr 



m=2 



and /C m (r 1; Ti; r m , r m ) is the m-th cumulant correlators, defined recursively 
by 

/Ci (ri,ri) = (v?(ri,r 1 )) / , 
/C 2 (ri, n; r 2 , r 2 ) = (T r [y? (ri,n) <p (r 2 ,r 2 )]) ; - /Ci (ri, n) /Ci (r 2 , r 2 ) , ... (6) 

etc. Further we will consider only the cases where /Ci = 0. 

To estimate T and electron energy 8 = lim^oo JF we use the same trial action 
as was proposed by Feynman for the polaron problem [24] S t = So+S pot where 

Spot = %f f [r (r) - r (a)] 2 e'^drda, (7) 
2 Jo Jo 

C and w being trial parameters. Then the Peierls-Feynman-Bogoliubov in- 
equality reads 



1 

where 



J~ <J~'t + — (Si n t ~ S pot ) t (8) 



T t = -1 In / e~ 5 ^ [r (r)] , („4> t = / „4 [r (r)] [r (r)] , 

/? Jr{0)=r(fi) ir(0)=r(/3) 

(9) 



which is equivalent to 



2/3 



■P r/3 m 

, , ... / (/C m (r (n) , n; r (r m ) , r m )) t J] dr 3 

m=2 m -P J0 J0 j=l 



oo m 



(10) 



To proceed, we will pass to the Fourier transforms 
(K, m (r (n) , n; r (r m ) , r m )) t 

= J ■■■ J Yl fcm(K 1 ,iuj 1 ;...;K m _ 1 ,iu} m - 1 )exp 



m— 1 



2 w i ( r i ~ r ™) 



x ( exp 



m— 1 



"X 1 9, D d D Kj 



(11) 



where Kj are the wave- vectors, VL D is the unit lattice cell volume, and a;j are 
the bosonic Matsubara frequencies. For the Gaussian trial action S t one has 



m— 1 



(^exp K i • [r (Tj)-r(Tm)} 



where 



exp 



1 



m—l 



Y f (Tj - T m ,T k - T m )Kj -K k 



j,k=l 



(12) 



/ (tj - r m , r k - T m ) = — ([r fa) - r (r m )] • [r fa) - r (r m )]) t = 

^ {([r fa) - r (r m )] 2 ) t + ([r fa) - r (r m )] 2 ) ( - ([r fa) - r fa)] 2 )J (13) 

Substituting Eqs.(17), (11)-(13) into Eq.(10) we find an exact upper-bound 
estimation for the free energy as a series in the coupling constant 



F < Ft — -p (Spot) — 

oo m ( . r ff 

£ " 

ra=2 



(5 m m\ J Jo 



^2 /C m (Ki,ia;i; ..;K m _i,ia; m _i) 



x 



m— 1 



o;i...u m _i 
-i m—l 



exp < Z ^ fa - T m ) - - J2 f ( T j ~ T m, Tk ~ T m ) Kj ■ K k } 
{ 3=1 3,k=l 

m ~ l n D d D K 1 dr j , 



(14) 
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In this paper we restrict ourselves to Gaussian approximation, which will mean 
ad hoc the neglect of the cumulant terms with m > 2 in the series of Eq.(14). 
Unless ip (r,r) is a Gaussian field indeed, the Gaussian approximation is be- 
lieved valid in a range of small enough g, necessarily satisfying the condition 



where W is a measure of the electron band width. Explicit criterion for ap- 
plicability of the Gaussian approximation depends crucially on the critical 
exponents and space dimensionality, see Section 3. 



3 Quantum case 

It was demonstrated by Feynman [24] that at (3 — > oo 

1 <[r (r) - r (,)]»), = ^ (l - e""'™') + ?£|r-,|,.W + ^ 
jj \ 1 1 v J v z w 

(16) 

(D is the space dimension) and so, with the notation A = v/w, we obtain 
and 

/ (tj - r m , r k - r m ) = (l - e -^-^l - e -«l^-r m | + e -h- fc |) 

A 2 

+ y (l r i ~ r ™l + l r fc - T m\ ~ \Tj ~ r k \) . (18) 
8 



Using Eqs.(16),(17),(18) and the Debye approximation for integration over K, 
to obtain 

Dv(l - X) 2 n r K ™»* ^ 
T < UVK \ A) - g Ad / £ /C 2 (if, iu;) x 



jf'cosur (l - j\ e-^-^-^) K2 dr 



K dK, (19) 



where A D = rjri — \ ' ^ ( x ) being the gamma function and if max is the 

Debye wave-number cutoff satisfying AuK® ax = D. For completing the limit 
P — > oo in Eq.(19) we use the method of residues to sum over the Bosc 
frequencies and employ the spectral representation 

K 2 (K, iu) = - I" J ( K,u) du, (20) 

7T J-oo U — ILJ 

with J (K, x) being an appropriate spectral density. So we obtain the varia- 
tional upper-bound estimation £ < £ (v, A), where 



£ (v, A) 



Dv (1 - A) 



2 



/ / / J"_(if,u)e T +2 J 2 » ^ ] K D ~ x dKdudt 
2n Jo Jo Jo 

(21) 

and J7"- (if, u) — J (K, u) — J (K, —u). Note that for even-frequency spectrum 
fluctuations (in particular static ones) J71 (K,u) = 0, so the interaction term 
in Eq.(21) vanishes. 

3.1 The use of scaling 

Until now the statistical properties of the field (p (r, r) have not been specified. 
Further we will use the dynamical scaling law near the QCP [2] 

X (K, u) = f^J_ (fK, f z u) , V/ > (22) 
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where r] and z is an "anomalous-dimension" and dynamical critical exponent, 
respectively. Using in Eq.(22) / = K ms , x /K, we have 

j -<*»H£r^M¥)'»)- (23) 

Plugging Eq.(23) into Eq.(21) and using the substitutions for the integration 
variables 



u 



K 

zu, K = K max ^, t = v^s, (24) 



K 

notations for the parameters 

W = ^K 2 nmx ,q=^,d = D-2 + V , (25) 
W being just the band width in the Debye approximation, and for the function 



( s , A 2 ) = (l - A 2 ) (l - e- s ) + X 2 s, (26) 
as well as rescaling g to fix the normalization of the fluctuation spectrum 

roc I — 

/ Q (zu)dzu = 1, Q (w) = -J(K max ,zj) , 

J 7T 



we obtain 



' 1 " ' ' * x 2 e L 



dxdsj 
(27) 

where the indexed by to angular brackets mean averaging with the weight 
Q{w). 

3.2 Weak-coupling regime 

In the weak-coupling regime /i — 1 — A 1, while the range of the parameter 
q is not predetermined yet (however, the restriction q < 1 should be imposed 
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anyway, otherwise the continuum description could not be used). In this regime 
the electron is weakly "fluctuation-dressed" . Using assumed smallness of \x we 
can expand the right-hand side of Eq.(27) in the Taylor series with respect to 
\i. This gives up to the terms of second order in \i inclusive 



p ( \\ D (A \9* D (A ^ 2 D 
S [v, A) ~ -— a (d, z) yy~j a i (d, z, q) j 



g,2 



q+ [yy] a 2 (d,z,q) 



Wfi 2 
(28) 



where 



/ 1 d+z i \ 
a (d,z) = I [ X \ L dx) , (29) 

WO X + §7X2 / 



W i vj 



d+z 



ai (d, z,q)=q( I — dx ) (30) 

\ Jo (x + §X2) (q + x + §X2) I 



and 



I n (2g + 3x + 3fxf)x^ 

a 2 (d, z, q)=4q ( - T - — 

V° (x + §x- 2 ) (q + x + §X2) (2q + x + §X2) 



dx 



-ai(d,z,q) (31) 



The first term in Eq.(28) is the electron band edge shift in the lowest-order 
Born approximation, the second term is the potential energy and the third 
term is the renormalized kinetic energy. 

The Eq.(28) is to be minimized with respect to fi and q. Let the optimum 
values of the variational parameters be /x and q . Within the small fi regime, 
the correction oc g 2 to the bare kinetic energy that describes the fluctuation- 
driven renormalization of the electron effective mass, results in a contribution 
oc g & to the optimal bound S . This contribution is negligible when expanding 
So up to terms oc g 4 inclusive. The condition that allows to neglect the above 
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renormalization reads 

g \ 2 \a 2 (d,z,q ) 



« 1, (32) 



which is, in general, consistent with Eq.(15). Assuming the condition of Eq.(32) 
to hold, we minimize Eq.(28) first in \x and next in q. This gives the following 
expression for /i and £ 

^ ~ lw) q (33) 



and 



S = -ja (d,z)^ - j ai (d,z)^, (34) 



respectively, where the positive number a± (d, z) is the maximum of the func- 
tion 

GWz ,,) = !Mm), (35) 

viz. 

di (d, z) = max G (d, z,q) — G (d, z, q ) , (36) 

0<g<oo 

and q is the point where this maximum is attained. Note that lim^oo G (d, z, q) 
due to Eqs.(30), (35), so for existence of the above maximum it would be 
sufficient that lim 9 ^ G (d, z, q) = 0. As deduced from the very structure of 
S (v, A) (Eq.(28)), the parameter 

1 (w\ 2 i 

h = V -= = — [ai (d, z)]~* K m l x (37) 

is a measure of the fluctuon potential-well size, which should be much larger 
than the lattice constant, i.e. satisfy loK max ^> 1. 

By the virtue of Eq.(31) |a 2 (d, z, q)\ > a\{d,z,q). Therefore, once Eq.(32) 
is checked to hold, it automatically results in /j <^ 1, due to Eq.(33). On 
the other hand, inability to satisfy Eq.(32) would mean inapplicability of the 
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perturbational regime. After this general analysis, let us consider different 
cases regarding the critical exponent z. 

3.2.1 The cases with z > 2 

In this case we always have <C x, due to smallness of non-adiabaticity 

parameter so Eqs.(29) - (31) reduce to the functions of the combined index 
d* = d + z - 2 

a Q (d,z)~A (d*) = ^, (38) 
ai (d, z, q) ~ Ai(tf, g) = g^$ d * (g) , (39) 

where 

-i -—i 

®b (x) = I" p-rrdt, b > (40) 
Jo t + 1 

and 

a 2 (d,z,g)~,4 2 (cr,g) = (11 + 2d*) A x (d*, q) - 8A 1 (d*, 2q) - j-j—. (41) 

l + q 

The necessary condition for the finiteness of the above integrals is d* > 0. One 
can see that in this case the fluctuation spectral density shape is completely 
irrelevant. 

To infer on existence of the maximum of G (d, z, q) = G d * (q) = g '*" 1 [$ d * (q)} 2 
we first note that lim g _^ (q) — at 2 > d* > 1, since for such d* 

lim$ d * (q) = r(-)r(l--)= (42) 

For d* = 2, $ 2 (g) = In (g- 1 + 1), and lim^ G 2 (g) = lim^ g In 2 (g- 1 + 1) = 
also. The functions $d* (g) for (rather unrealistic) case 4 > d* > 2 are 
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reduced to those with d* < 2 using the functional relation 

, d* > 2 (43) 



d* 

q-$ d * (q) = q 



2 d*-2 

5 2 ®d*-2{q) 



d* -2 



and again we get lim g ^ Gd* (g) = 0. As outlined in the previous subsection, 
this means that at least one maximum point < go < oo does exist at d* > 1. 
One the other hand, the equation J^G^* (q) = for determining g is rigorously 
transformed to the following one: 

2-d* 

(<T - 1) = 0, (44) 

which obviously has no solution if d* < 1 . Thus for < 1, the weak coupling 
regime never applies. This exponents range will be revisited in Section 3.2. 

For 1 < d* < 2, the assumption of small g would allow one, by the virture of 
Eq.(42), to solve approximately Eq.(44) in a closed form. However, compared 
with numerics for specific d*, this approximation seems to be too inaccurate. 
An approximate equation, which results from inclusion of the next-to-leading 
terms of that asymptotic, can not be solved analytically anymore. So given 
d* , a reliable calculation of go requires numerical approach. For some cases of 
rational d*, one of them is considered below, $ d » (g) is expressed in elementary 
functions. 

Let us put d* — §. This case is a representative for fractional-rational d*. We 
have 

M = In — y + \ h arctan (V2q~* + l) + arctan f \plq 

V2 J q -2 + + 1 v ' v 



G. (g) = 2g5 



In — — \ 1- arctan f v^g * + l) + arctan f v^g * 

g"2 + V2q~* + 1 V 7 V 



The graph of G| (g) is shown in Fig.l. Eq.(44) for d* — | has unique solution 
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Fig. 1. Graph of the function G3 (q) 

2 

go ^ 0.126, for which Gs (q ) = ai (§) ~ 1.589. Checking Eq.(32) yields after 



cumbersome calculations 



W 

Provided that Eq.(45) holds, we obtain from Eq.(34) 



^7 < 0.378 , (45) 



and from Eq.(37) 



/O^max ^ 0.793 Ojj . (47) 



The numerical results obtained for different values of d* > 1 show that within 
the weak copling regime the smaller d* the larger numerical factor of the 
fourth-order correction in So, and the narrower the range of g where that 
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approximation works. 



3.2.2 The cases with < z < 2 



For < z < 2, Eqs.(29) - (31) are transformed quite specifically. Let vj$ scales 
fluctuation frequencies, so that Q (zu) be a function of the reduced frequency 
v — Then, omitting from now on the index of the averaging over vo (or 



u), we have 



U / 2-z \ — V W 



(48) 



and $ & (x) is defined by Eq.(40). It is seen that in the present case the weak 
coupling regime has a sense only at d > 0. 

The asymptotic of a (d, z) at ^ <C 1 depends critically upon the sign of 
d + z — 2, yielding 

r / a 

1-7 



do (d, z) ~ < 



, -, 3 > W a -, d+z-2<0 



(49) 



d+*-2>0 

where lnz7 = (lnz/) and Eq.(42) is taken into account. Thus, the Born energy 
scale depends on the fluctuation dynamics: (i) drastically in the first subcase, 
including in particular original Feynman's polaron [24]; (ii) weakly in the 
second subcase; (iii) negligibly in the last subcase. 

Next two integrals (30), (31) are transformed and asymptotically represented 
at — 3> 1 as follows 

G70 



ai(d,z,q)~[ — At (d,z, x) , i = 1,2, 



(50) 
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2 

where x = q 2 ~ z is a new variable to optimize over, and 



A 1 (d,z, 



2x 



u 2 -zdu 



2-z\Jq (u + v) 2 (x + u*-* + vu*-^ I 



(51) 



A 2 (tZ, z, x) 



SX" 



2 - z \ Jo 



u 



- A 1 (d, z, x) 



(2x + 3u 2 -* + 3z/w 2 - z ) u 2 -* +1 du \ 
+ u) 3 {x + M~ -|- v U ~^j (2x + -u~ + vu~^j I 

(52) 



Note that for any reasonable d the integrands in Eqs.(51) and (52) fall off 
at u — > oo faster than u~ 2 . Therefore, in the both integrals, unlike that in 
Eq.(48), the upper limit ^ ^> 1 has been safely replaced by oo. 

In the case considered the expression (35) is parametrized as follows 



G[d. .:.(/) = (— V 2 ~* G(d,z,x), G(d,z,x) = A i^ z ^\ (53) 

^VJqJ X 



Accordingly, the energy asymptotic in weak coupling regime at z < 2 is given 
by 



AW \wj 



for z < 2 - d, 



2 {V 



g \ 2 



2-sin(g) U» ^ } 



(54) 



for z = 2 — d, and 
for z > 2 — d, where 



2, / W 
— In 



d + z-2 \W J W 



_ o 3 — 2 — d 



(55) 



(56) 



(d, z) = max G (d, z, x) 

0<x<oo 



(57) 
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Finally, Eq.(37) yields for the fluctuon size in the present case 

l,K^^(^j(^f^ Wd,*)]-* . (58) 

For a given g, the fluctuon size at z < 2 proves parametrically much smaller 
than that at z > 2 unless d + z > 3. 

The sufficient condition for the perturbational regime to hold is provided by 
Eq.(32), which reads in the present case 

\A 2 (d,z,K )\ ( g\*(W\" 



2-z 



< 1, (59) 



x \W J \zu , 

where the numerical factor requires a numerical calculation for specific d and 
z. This condition proves much more stringent than Eq.(15), but assures that 
I i^ max ^ 1 in any case. 

Now the key question is that of existing the optimal x , to answer which 
exploring the behavior of A\ (d, z, x) at x — > is crucial. Let us assume that 
(u- a ) < oo for all a > 0. Then at z < 1 + ±d 

Ai(ci,z,x) 2-d / _4-d-«\ (^frf) 

hm = 7t(is ) — — -At- < oo, 

x (2 - \ / S in (vrfef ) 

so lim^o C (d, z, x) = 0. At z — 1 + \d we obtain the asymptotic at x — > 
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3 2— d 



2 + a x 7 x / 

so we have lim^_>o G (d, z, x) = also in this case. Hence at z < 1 + \d the 
maximum point x surely exists. At z > l + \d, making use of the replacement 

2—2 2—2 

u = (xz/ ) z £~ and of Eq.(42), we obtain the following asymptotic 

2 d+2-2 / d+2+z \ 7T 

Ai(d,2, x ~ -x * (z/ - ) / ttz \ i x — » 0, 

2 \ / gin ^Tr^+I^j 
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from which we infer that x exists, since lim^ G (d, z, x) — 0, if z < | (d + 2) 
(that holds automatically for d > 1). If z > § (d + 2), which may occur for 
< d < 1, the above limit is either a finite number or oo that makes weak 
coupling regime nonexistent. 

For completeness, it is instructive to consider numerical examples. We consider 
two important cases z = and z — 1 falling into the class z < 1 + \d, for 
which the existence of x has been proved above. In the both cases the relevant 
formulas, before the v averaging, are expessed in elementary functions. Due to 
persisting v averaging and arbitrary d, however, the formulas yet remain too 
comlex for illustrative numerics. To make things simpler, in the subsequent 
two examples we assume that d — 1 and the v distribution is strongly peaked 
at v — 1. We do not expose the corresponding graphs of G(l,z, x) since 
they are pretty much similar in shape to the graph shown in Fig.l, apart of 
appreciable difference in scales of variables x and q. 



Example: d — 1, z — 0. With the above assumption this is actually the 
Feynman polaron problem [24,25] . We obtain 

G^M^iW-^^) 2 - (60) 

X \Z H J 

This function achieves its maximum at x = 3 in accordance with Feynman, 
which gives A 1 (1,0) = f^. Then Eq.(54) reproduces the Feynman result for 
the energy bound 



6 = ^„4, a ( 61 ) 



81 r 4 W / \w 
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while Eq.(58) yields the fluctuon (polaron) size parameter in terms of Feyn- 
man's a constant 

k>* _ ^ - • (62) 



6\/6 y/m^ V81 

These results have a sense upon satisfaction of Eq.(59), which now reads 



(63) 



Example: d = 1, z — 1 This case corresponds to the interaction with 
acoustic-like critical mode. Now, one should maximize the function 

/-. -i x 4 /3x — 1 arctan — 1 x— 1 \ 2 

G(l,l,x) = - ; + lnx-1 . (64) 

V ; x V x V4x- 1 2x / v ; 

We find numerically that the unique maximum point is x ~ 3.81 and A\ (1,1) - 
0.208. Then, Eqs.(55) and (58) yield 



2 



c 3 9 

Cn ~ 

° 21f 



In ( — ) +0.104 



, ...-?-) 



(65) 



and 



/O^max ^ 2. 19^^, (66) 

9 

respectively. In the present case the condition for the perturbational regime, 
which doesn't contain W at all, reads 

i^(M.*)i m' g0 . 112 (67) 

xo Vtuo/ Vtuo/ 

or \g\ <C 3zuo. 

To conclude this section, for < z < 2 weak coupling regime is realized at 
much smaller g than for z > 2. For the latter, g should fit Eq.(15) while the 
characteristic fluctuation frequency ro plays no role. For the former, however, 
the upper bound of \g\ /zuq, is crucial. 
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3.3 Strong coupling regime 



3.3.1 General Consideration 

In strong-coupling regime the electron is heavily "fluctuation-dressed" . Let us 
make in Eq.(27) the variables replacements y = g _1 x, r = ye~ s and t — 1 — r. 
This transforms that equation to the following one 



S (q, A) = ^Wq (1 - A) 2 - j^- 1 (M (q, w, X- 



where 



(68) 



M (q,w, A : 



d+z 



y 2 e 



dy 



« /'a^i -(i-«> 



+ ( 1 - A2 )i *i 



o e(y,g,ro,A 2 ) 

e^^y^y, 



and 



/o e(y,g,ro,A 2 ) 



e [y,q,w,X J = — g^ y? + X y. 



(69) 



(70) 



To proceed, it is important to note that the function e (y, g, zu, A 2 ) increases, 
in the integration range over y, from zero to g _1 e, where e = ^ + A 2 . Thus, at 
g ^> e Eq.(69) may be expanded in asymptotic Laurent series in overall small 

e(y,q,w, A 2 ) 

oo 

M(g,^,A 2 ) = £ M p (g,zi7,A 2 ), (71) 

where 



P =-i 



M p (g, n7, A 2 ) = 
with A^_! (0 = e"*, 

^(0 



6 (y, g, «7, A 2 )] P iV p [(l - A 2 ) y] y^dy (72) 



-if r l 



(p + iy.Jo 



e~* [f(t)] p+1 t^dt,p>0, 



(73) 
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and 

l n (i — f) 00 f k 
/(t ) = _Ml_fl = Ss L_,0< t <l. (74) 

The Taylor series representing / (t) converges at [0, 1) and so does the Taylor 
series for any integer power of / (t) 

00 

[f (t)Y l = E a n,mt m . (75) 
m=0 

Typically, the strong coupling-regime fTuctuon binding energy ~ Wq is smaller 
than the fluctuation energy. Hence the above-assumed relation between q and 
e is satisfied if A 2 <C q. Another, weaker, criterion for expanding M p (q, zu, A 2 ) 
in powers of e (y, q, zu, A 2 ) is inferred on by noting that a left vicinity of t — 1 is 
the dominant range for the integration over t in Eq.(69). Hence at A 2 <C 1, it is 
the range y < 1 that contributes mostly to the corresponding integral over y. 
In this range e (y, q, zu, A 2 ) < ^q^' 1 + A 2 , is small, unconditionally for z > 2, 
and under the condition Wq 1 ^^ ^> zu for z < 2. Actually, when truncating 
the series of Eq.(71), either A 2 q or A 2 1 and Wq l ~^ 3> zu are our the 
only approximations. We should check them at the end of our calculations. 

Let us try to simplify the above-developed expansion, by picking in it up the 
leading terms with respect to k = (1 - A 2 ) not imposing in advance 

any other restriction on q and A 2 . To this end let us transform M p (q, zu, A 2 ) 
as follows. For M_i (q,zu, A 2 ), we obtain directly 

q—- 1 M_ 1 (q,zu,\ 2 )=Ki —du, 76 

Further, using in Eq.(72) the Newton's binom, we arrive at the identical but 
more convenient representation 

q^M p (q,zu,\>) = (l - {^f A 2 ^>m p+1 , fe+1 («) (77) 
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where 



m nj i (k) 



(-1) 



n— 1 



/ 

JO 



1 7 (id, + n + 2,«) 



[/(/™)] n dw, d, = d+ (z-2)l 




(78) 



(di = d* which has been introduced in Section 2 for the case of z > 2) and 



is the incomplete gamma-function. The integral in Eq.(76) at z < 2 converges 
if d > irrespective of A, while at z > 2 this is so if A = strictly. For A ^ 0, 
even small, the convergence condition at z > 2 reads di > 0. These restictions 
upon the critical indexes are the same as in the weak coupling regime. 

The value of M_ x (q, zu, 0) is independent of z, and given by 



However, estimating M_ x (q, zu, A 2 ) at A 2 ^ 0, except for the case z = 2 
where the factor (|^ + A 2 ) plainly replaces ^, depends crucially upon z. 
We postpone this task to consideration of specific cases. At the same time, 
asymptotic series in k for M p (q, zu, A 2 ) with p > can be obtained by an 
independent of z trick. 

Substituting the series of Eq.(75) into Eq.(78), integrating by parts and using 
the well-known asymptotic 




(79) 




(80) 



7(6,x) = r(6)+0( ; 



x 



x 



) 



, X > 1, 
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we obtain with an exponential accuracy 



m»,i («) 



: _ ir -i 



{ 1 \ 



'm + n)la. 



n,m m 



+ c n>l (m t + n)\ fin k 1 - (m t + n) — ) n mi 

V mi + nj 



(81) 



mi + n, 

where c n> ; = a njm; , mi being an integer, if any, satisfying the condition 2m ; = 
di, and otherwise c ni 2 = 0, 

bn,i = t ~^h-, (82) 
j m — Mi 

and ^(x) is the logarithmic derivative of the gamma-function. At < d < 2 no 
mi emerges if z = 2, the only m± — may appear, if z < 2 (e.g. for z = d = 1), 
and if z > 2 an infininte number of mi > 1 may exist for some d. It is worth 
noting that at z < 2 and di ^ 

Kt = f \mr_i dt _ i. (83 ) 

Only the cases with mi — 0,1 may be important since the O (n m In terms 
with m > 1 are small compared to the kinetic-energy term in Eq.(68). By the 
same reason, of the series in integer powers of k in Eqs.(81) we retain only 
the 0(1) term that exists unless mi = 0. The thus approximated Eq.(68), 
after performing some interim summations over p and neglecting purely non- 



adiabatic corrections O y j, becomes 

So (q, A) ~Wq (1 - A) 2 - ^g^¥ (w^U A 2 ) 



4 W 

where 



A(q,\ 2 )+A, (84) 



A = -£iw {l - s '"° ] (85) 



24 



is an energy shift, independent of the variational parameters, 




«2 e 



o 1 + X 2 x~ 1 u 1 ^ 



oo n I I \n— 1 /I \ _ 

E E ^-in— <^ir (\di + n + l) b n:l (l - A 2 ) " X^x 1 (86) 

TI=1Z = 1 VZ 7 



and 




In (1 - A 2 

X 2 



+ 5 du2 (l-A 2 )g(ln-+ 7 -- 




with 7 being the Euler constant. In Eq.(87), the first and the second term do 
not emerge at z > 2 (where di > necessarily) and at z < 2, respectively. 
In all cases where d\ > 0, A = — £ B , the band-edge shift in the lowest-order 
Born approximation. 

Further analysis on the base of Eqs.(84) - (87) depends crucially on whether 
z > 2 or z < 2. We consider these cases separately, detaching z — 2. The 
peculiarity of the latter case allows us to calculate II (x, A 2 ) in a closed form 
and, that is not feasible in other cases, to ultimately explore an impact of the 
spectral weight Q (zu) on the fluctuon formation. 

3.3.2 The cases with z = 2 

For z = 2, di = d and b n> i = b n>n , so that A (q, A 2 ) = and Eq.(86) greatly 
simplifies. The answer reads 



£ (q, A) ~ (1 - A) 2 q - ^ 9 -Y Qd) (l - A 2 ) i?, (A 2 ) - (88) 



where 




(89) 
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and h{t) = f(t) — 1. Formally Eq.(88) matches the case of A = 1, as the 
fluctuon binding energy obtained vanishes at A = 1, that concords with exact 
Eq.(69). However, that point is likely isolated since in essential weak-coupling 
regime, i.e. at < 1 — A <C 1, the condition at k <C 1 may break down. 



Minimization Eq.(88) first in q and next in A, we find the optimal q value 



'd 



I ( t - 1 I ( R d (Xl 



WJ 



1 + Aq 
1-Ao, 



(90) 



as well as the bound energy 



e D ( 2 i 



WJ 



2 

2-d 



W-S B , (91) 



where 



P d (\) = (l + \)(l-\) 1 - d R d (\ 2 ) (92) 
and Ao is the maximum point of the function P d (X). For d ^ 1, Eq.(91) 
presents a singular perturbation expansion in coupling constant. When A 
corresponds to an extremum, it satisfies the equation 

R> d (X 2 ) d-(2-d)X 

2X MX^) + i-a* " °' (93) 

otherwise A = 0. For the latter case, 

P d (X )=R d (Xl) = — + 0(1), 

where w = (zu^ 1 ) \ which attains, to within 0(1) terms, largest of all 
possible values of those functions. Note that \im m ^ Q (zu) = 0, so it is likely 
that < oo. 

Let us search a solution A to Eq.(93), in the vicinity of A = 0. Assuming that 
also (w~ 2 ) < oo, we have in the leading approximation 



R d (X 2 ) ~W (w- 1 ) - X 2 W 2 (w- 



(94) 
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which yields for the sought solution 

A » * 2W "* = (95) 

For "rigid" i.e. zeroing below some finite the above-exploited as- 

sumption (ro~ 2 ) < oo holds automatically. Consider now "soft" Q(zu), for 
which (tu~ 2 ) = oo, but (zu^ 1 ^ a ) < oo with some < a < 1. Scaling the 
behavior of Q (-07) at w — > + by 

^ / n , / i „\ sin (7TO") _ , 

Q (tu) ~ 6 CT ) — -w CT , b a = const, 

^ ' 7TO" 

we obtain the solution to Eq.(93) at 1 > a > 1/2 

X ^{2k) (w) ' ro '=((^)j ' (96) 
If < a < 1/2, A remains zero. Since ^zr > 1 at 1 > er > 1/2, the non- 

a 

adiabatic corrections resulting from A ~ (yf^J 2(7-1 are even smaller than those 
~ ^ resulting from the integral term in Eq.(89). 

Thus, as far as small A is concerned, either Ag = o or A = for all 
admissible Q (w). Neglecting the postleading non-adiabatic corrections, from 
Eqs.(90) and (91) we arrive at 



2 I Ww 



2 

2-d 



' £ o = -j(l~ 1 ) Wqo ~ £B ' (97) 



Requiring q <^ 1, one gets the criterion of applicability of the continuum 
approximation 

r (H^ <L < 98 > 

Under this condition, the self-trapping term oc i^yy) 2 d in £o ma Y be both 
smaller and larger than Sb- The latter situation occurs if coupling is strong 
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enough to satisfy 



I 2 I Ww 



> 



2 2 1 



d2-d 



r(|d) 



2-d 

d , 2-d 
G7 



(99) 



Even though — £ B dominates £ , the self-trapping term yet lowers S more 
than does the correction oc (jys) in weak coupling regime. 

Consider now the singular case d = 2, for which m\ = \ and 
Here we obtain from Eq. 



= /M/Wl--i-^_ L 

JO 



t 2 



£o (9, A) = 



W(l-A) 2 -^(l-A 2 )i2 2 (A 2 ) 



+ 



where 



ifc (A 2 ) 















) = 










Jo 



h(t) 



+ 



h(t) 



dt 
t 2 , 



(100) 



.(l + e/i(i)) 2 ' l + eh(t) 
This expression is easily optimized first over q and afterwards over A to yield 



g = e 3-7+S(Ao) > ^ = _( 1 + go )^ ) 



(101) 



where 



-(7)1 



/W\ 1 - A 
+ A' 



and A is the maximum point of the function S (A). Searching again A <C 1, 
we obtain 

Ao^ 2/ 2X , go = e wo ^ (102) 
It is seen that for Ao 1 and qo <C 1, the inequality g 2 (ro~ 2 ) 3> 1 and Eq.(98) 
with d = 2, respectively, should hold. £ given by Eqs.(lOl), (102) is much 
above that obtained in weak-coupling regime (see case d* = 2 in previous 
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subsection) for typically -^p -1 = (1), though if ^ -1 = (sa) the 
former may gain. 

But what happens if Eq.(98) doesn't hold ? The answer is easy for 2 > d > 1 - 
in this case weak coupling regime may realize. For d < 1, however, the question 
cannot be answered within the present framework, as numerical study reveals 
no any maximum of Pd (A) other than that in a close vicinity of A = 0. 



3.3.3 The cases with z ^ 2 

Using the experience with z = 2, in what follows we restrict ourselves to small 
A, and assume (w~ l ~ s ) < oo, where < s < 1 throughout. The integral part 
of II (x, A 2 ) possesses small-A expansion at A 2 <C x, which we force to hold. 
Further, we have x C 1 unconditionally if z > 2. For z < 2 we force holding 
x <C 1 anymore. At the end, we check those conditions both to hold. With 
such prerequisites, up to the first-order terms inclusive, we obtain 

S (q, A) = S (q, 0) - ^WqX + ^WT (l + (^ 2 ) g 1+ ^A 2 , (103) 

at d > z — 2 and 

e (4 , a) = £o M - > 9 a + ^^>1Z: (KY \ 2 " (m 

2 2d sin ira p (l + |) W CT / 

at d > 2 — 2, where <r = Here 



£b («, o) =£ W , - g [9 (0)] V ?i - _ I [? (0)] V ^^y^" 

+ A-^A(g,0), (105) 
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q (0) is g obtained with A = 0, i.e. given by Eq.(97) and 



6 = 61 

x 1 



z + d 



(l - sM) + 2 7 , z + d + 4 
1, 2 + d = 4 



(106) 



Let d\ 7^ 0, 1, i.e. A (g, 0) = 0. For d > z — 2, the minimization equation for A 
is solved to give 

z — d 

q 2 



A(?) 



(107) 



s 2 (sr- 2 >r(i + ^)' 

Then the minimization equation for g is well solved by iterations in small 
adiabatic parameter, to yield for the variational parameters 



where ac? /10c is defined by 



(1 - *±*) r & _ 1) r (1 + f) 



r(i + f) 



r(i + ^) 



(110) 



For d < 2; — 2 that may realize only at z > 2, we find the optimal A at a given 
q to equal 

1 

(^J [g(0)] — — g— , (111) 

which provides a minimum if 2<r > 1 (i.e. d > £ ^), otherwise we should put 
A = 0. Just as above, the equation for optimum g at the present conditions is 
solved by iterations, which results in 



\ / \ z — 2 singer / d 
AW = — — r 1 + 2. 
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9 = 3(0) -^^[„(0)]» (112) 

^(^MO)!*-)*. (I") 

where here w denotes only the first term in the expression given by Eq.(llO). 
To check all necessary conditions, we consider below the cases with z > 2 and 
z < 2 separately. 



z — 2 sin tv a ^ ( d 
A ~ r 1 + - 

2 7TO~ I 2 



Subcase z > 2 For z > 2, we have from Eqs.(108)-(113) g = q (0) + o (§) 
and A = o Both A n 1 and Aq <C ^q§ are satisfied automatically. 
So the corrections to formula for £ as given above for the cases with z = 2 are 
much smaller than w and even not worth to be considered anymore. There 
remain the same conditions, given by Eqs.(98) and (99), as with z — 2. 



Subcase z < 2: For z < 2 and d + z — 2 7^ 0, using Eq.(108) we have for 
original parameter v = qW up to the first order corrections 



v = q W ~ W 



9' 



v 2 + I Wtu 



2-d 



where, as introduced above, 



w = 2 



(i-^)r(^)6 (^-i)r(i + f 



r(i + f) 



w + 



r (1 + 



^1, 



for the parameter A 



A n 



r(i 



2 



r(i + f) 

and for the fluctuon energy 



\2 1 Wzu 



2-z 
2-d 



£o = ~(- 1 -i)w 



4 \d 



I 2 I Wzu 



2-d D 

+A -4 



w, 
(114) 

(115) 



w 



\2 J Wzu 



(116) 



;ii7) 
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where 

r(£&)b r (1 + 1) 

w oi = —7 — h\ w ® H — r — r^^i ( 118 ) 
r(i + f) r(i + ^) 

Now the check of necessary conditions is in order. If we require that O (X 2 x~ 1 ) 
terms should be small on average, we get 



r i 



d-z 



2-d 
2-z 



(119) 



r(i + |) 

The conditions that i C 1 and A <C 1, to within purely numerical factor, 
give the same inequality as Eq.(119). 

Note that at d+z—2 < the value A > and has no connection to £ B . In these 
subcases, Eq.(119) proves much stronger than that of Eq.(99) that leads to 
total domination of the self-trapping energy term over A. Moreover to within 
the present approximation, A is much smaller even than the O (zu) correction 
in Sq. As an example of such a case, consider again Feynman polaron (D = 3, 
d = 1, z = 0). From Eq.(106) we have b = 41n2 and from Eqs.(114) - (116) 
we obtain, in terms of Feynman's a, for the original variational parameters 
v = Wq and w — Xv 

(Aa 2 , \ 

v ~ h 1 — 8 In 2 w , w ~ zu 

V 9tt / 



and for the energy 



^ = -g + 61n2 + ^ w , 



These results are valid upon the conditions 



— > a > ^V27rln2 ~ 3. (120) 
2 y zuq 2 

The left-hand side inequality (particular case of Eq.(98)) doesn't appear in 
Feynman theory, since there W = oo. 
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At the end, explore singular cases, with mi = 0, i.e. d + z = 2. Using Eqs.(84) 
- (87) we obtain 



£o (q, 0) = jWq - jWq\ + ^WT (d) g 2 (^ 2 ) q d \ 2 - 
-\ q m^qU- 9 -\n (S- 



(121) 



As above, the minimization equation for A is solved exactly 



A(?) = 



A-d 



while that for q = q (0) y, being 



1 — dw\ 
1 - tt^t— [q 



T(d) W 

is well solved by iterations around y — 1, to yield 

^ 9(0) + { b(0)] H_j^__2_£ 9 (0) } r ( 1 + £ 

Eq.(123) is valid provided that 

cA 2 — d W\ 



r i + 



[? (o)p<i, 



(122) 



2; r (d) w 

which means a sort of strong-coupling conditions, considered above 



2 

"2-d 



(123) 



r i + - ) > 



2-d / d\ 

r (d) \ 2 J V w 



(124) 



2 y ro iu 

Then using Eq.(121), (122), and Eq.(123) we obtain in the leading approxi- 



mation 



r( d - + i )J?- 

\2 ) Ww„ 



2-d 2 
+ - 



r(i + j)] 



Y{d) 



g Y m 



— — , (125) 



An ^ 



r (f + 1) 

T(d) 



\2 Wvj q 



d 

' 2-d 



ZUl 

W 1 



(126) 
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and the energy 



£ - D ( 2 



1 W 



9' 



Dg 2 , 



Wwq 

2 ' 



9 



r d - + i 

I 2 I Ww 



2 

2-d 



2 

2-d 



3D g* [r(l + f); 



4 W 



r(d) 



CE7i 

(127) 



As an example of the peculiar case d + z = 2 one may consider z = d = 1. 
Assuming for simplicity tui = E7 we obtain 



7T 



90 



w 



) 2 , A 



Wo 

, 9 , 



and 



£() = 7T 

16 



^) +3 

wo 



£ + l n f v^ 7 ' 1 ^ 
_ W 2 W 11 V 2 W^ 0/ 



(128) 



Strong-coupling condition (124) in the present example simplifies to 



(-) 2 » 1 



It appears that this condition and weak-coupling condition given by Eq.(67) 
have wide overlap, within which Eq.(128) results in much lower £ than Eq.(65). 
Even the absolute value of logarithmic correction proves larger than that of 
Born shift -ffp In (^)- This means that the strong coupling solution is ener- 
getically more favorable in the overlap region. 



34 



4 The self-trapping and electron density of states at classical crit- 
ical point 

4-1 Variational estimation for the electron free energy 

Let us consider now the self-trapping of the electron at a classical critical 
point, CCP (or second-order phase transition) at finite temperature T c = 
(rigorously speaking, the transition can be considered as a classical one 
only assuming that it is not too close to QCP at zero temperature [2]). The 
Feynman variational approach has been applied to this problem by us earlier 
[21,22,23] (only for a particular case D = 3, 77 = 0) but here we reconsider 
this (for a generic situation) concentrating on some new points such as the 
behavior of the electron density of states (DOS) and detailed comparison with 
the quantum case treated above. 

We start with the same general expression given by Eq.(14). Typically for 
CCP one has h~wf3 c <C 1 due to well-known phenomenon of critical slow- 
ing down [31]. This is true provided that a typical wave vector of the order- 
parameter fluctuations is small in comparison with the reciprocal lattice vec- 
tor; in our case the typical wave vectors K* ~ l/l (where Zo is an opti- 
mal fluctuon size) should be much smaller than K max and therefore, indeed, 
hwf3 c ~ (K* / K max ) z < 1 so we can use for our estimations long- wavelength 
asymtotic of static order-parameter correlators. Due to irrelevance of the dy- 
namics one can put it the trial action (7) w = 0. We will also use the notation 
C = u 2 /2 where uj is the frequency of the trial oscillator; the fluctuon size is 
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/ = (h/2muj) 1 ^ 2 . We will be interested in the strong-coupling regime where 

(3 c hw > 1 (129) 

Then instead of Eq.(19) we will have for the Gaussian case the following 
estimation (cf. Ref. [21]) 

Duj 3d 1 An f K ™* ( K 2 \ 

n-T-^l MA')exp (--)*-> AK (130) 

where /C 2 (K) is the Fourier transform of the static order-parameter correlation 
function with a small- K expression 

lC 2 (K) = (^pj (131) 

A numerical factor factor in the above expression is absorbed into the coupling 
constant g. For the reasons which will be clear below we consider (5 in the 
partition function and, as a consequence, in Eq.(130), a running variable. 
Substituting Eq.(131) into Eq.(130) one promptly finds 

After minimization of the right-hand side of Eq.(132) we find for the optimal 
estimation of the electron free energy 



DW{2 - d) 
W,9) = W — 



r (i +1 \f>9 



= - BW [w) (133) 



2 W 



Similar to Ref. [21] one can show that this is an optimal estimation provided 
that 

(pW) d/2 < (f3g) 2 < (pW) d (134) 

where the left inequality gives the criterion of the strong coupling, or self- 
trapping, and the right one gives the criterion of applicability of the Gaussian 
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approximation. The latter is found from the consideration of the scaling prop- 
erties of higher-order cumulants in the expansion (14). For (/3gY < {(3W) a/z 
(weak coupling regime) the second-order Born approximation turns out to be 
optimal. For d — 1 these results coincide with that from Ref. [21]. 

Comparing the result (133) with the ground-state energy estimations for strong- 
coupling regime (97), (117) one can see that in the leading order these expres- 
sions differ just by a natural replacement of the temperature fi^ 1 for the clas- 
sical critical point by a typical fluctuation energy for the quantum case. How- 
ever, the physical meaning of these quantities is essentially different: whereas 
for the quantum case we have derived an estimation for the true boundary of 
the electron energy spectrum, for the classical one our result is connected with 
the fluctuation density of states tail which is not restricted (in the Gaussian 
approximation) from below. Further we will prove this important statement. 

4-2 Electron density of states tail: Laplace transformation 

The electron partition function (2) can be estimated, due to Eq.(133), as 

Z ~ exp (BW^^g^) , (135) 

At the same time it can be rigorously expressed as a Laplace transform of the 
electron DOS 

N(E) = {5(E-H)) f , (136) 

namely, 

roo 

Z= / N(E)e-? E dE (137) 
Jo 

We can use now Eqs.(135),(137) to find the asymtotic of the electron density of 
states (that is why it was important to consider j3 formally as an independent 
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variable). Using the saddle point method one can prove that at large enough 
negative E 



N (E) oc exp 



2 \4-d 



'(§)' r -s 



2-d/2 



(138) 



with a suitable choice of the energy scale E as 

7lD 



2 

4-d 



E = 



4 d_ 

g4-<iW 4 - d 



(139) 



(origin of a numerical factor in the definition (139) will become clear in the 
next Subsection). The saddle point method is applicable if the exponential 
in the above formula is large, which is connected with the left inequality in 
Eq.(134). Another restriction is obvious from the observation that the real edge 
of the spectrum for the Hamiltonian (1) without fluctuation dynamics equal 
to E min = — gmax\ip\. Therefore the asymptotic (138) makes sense only for 
\E\ <C \g\. Near the edge of the spectrum the "Gaussian" tail (138) transforms 
into the "Lifshitz" one. Analyzing the scaling properties of the higher-order 
cumulants one can demonstrate that at E — > E min + 



N (E) oc exp 



const 



(E - E ■ ) d/2 

J -'min) 

This result has been obtained in Ref.[23] for d — 1. 



(140) 



4-3 DOS tail: diagrammatic approach 

To better appreciate the above-mentioned approximations, it is instructive to 
reproduce the result (138) by another way basing on the diagram technique 
[18,32,33]. The average Green function of the electron describing by the Hamil- 
tonian (1) with the Gaussian random static field ip (r) is written in a closed 
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form 



G (^ P ) £-P 2 /2-£(£,P) 

£ (E, P) = g 2 Q D [ 7 (P - K, P, K; E) /C 2 (K) P - K)-^5y (141) 

where £ and 7 are the self-energy and three- leg vertex, correspondingly, K, P 
are, as before, D-dimensional wave vectors, and static correlation function is 
given by the expression (131). To find asymptotic of DOS for large enough 
negative energies one can use a method proposed first by Keldysh for doped 
semiconductors [34] (the same trick was used also for magnetic semiconductors 
near T c [35] and for electron topological transitions [36]). For large enough 
\E\ , E < one can neglect momentum dependence of both E and 7 since only 
the momentum transfer K — > is relevant for d < 2. Also, we can express 7 
in terms of X via the Ward identity [32] 



7(P,P;0;£) = l-^|i3. (142) 



Then, taking into account Eq.(131), we obtain a closed differential equation 
for the self-energy of the form 

^ {E) -y--dE-) 9AD h E-K*/2-XW (M3) 



Consider now the density of states (DOS) 

An f K ™* K D ~ l dK 

Nd (E) = -— Im / — s ^— (144) 

y 1 7T Jo E-Z(K,E + i5)-±K 2 v ; 

It is clear that at \E — £ (k, E + iS)\ <C ^K^ iax at least for L> < 3 the main 
contribution to Np (E) comes from small K (K <C K ma , x ) region. 



39 



Let us solve now the equation (143). Integrating over K one derives 



E(E) = 



ttD g 2 

2 sin ^ w 



dZ (E) 
dE 



[Y,{E)-E]-> 



-l 



(145) 



./ 



E 

En 



(146) 



Denoting 

£ (£) — E = En 

with i?o given by Eq.(139) we obtain a non-linear first-order ordinary differ- 
ential equation 

For d = 1 this is Riccatti equation, which was solved in a similar context 
earlier [35]. We consider here only the asymptotic behavior of the solution at 
E < and \E\ 3> E directly from the initial equation (143). For these E 



|Im£ (E)\ < |ReS (£)| < |£| 



(148) 



and we linearize this equation with respect to the imaginary part of the self- 
energy to obtain 

dhnE(E) 1 f-E\ Y - 



dE 



Eq \ Ey 



ImS (E) 



(149) 



Thus, we have 



ImE (E) ~ CE exp 



4-d 



-En 



2-S 



(150) 



where C is an undetermined integration constant. At these energies, the den- 
sity of states becomes 



N D (E) ~ 



CD{2 - £>) £ 



2sin^ W f\ E 



2. i 7-i 1 2- ^ 



• cxp 



4-d 



E 



En 



(151) 



which coincides with the result (138), with an accuracy of a numerical factor 
of order of 1 in the exponent. This may be considered as a justification of our 
treatment basing on the Feynman variational approach. 
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The physical meaning of the self-trapping energy for quantum and classical 
fluctuons are essentially different. For the fluctuon near QCP, as well as for 
the Feynman polaron, we calculate approximately the ground state electron 
energy, or the edge of the spectrum. If we will calculate next-order corrections 
to the electron free energy in T = f3~ l we will find just a temperature shift 
of this energy rather than any exponential tail of DOS. The energy of the 
classical fluctuon is just a position of the chemical potential at small enough 
electron concentration n. For 



which is a capacity of the tail, the chemical potential level is "pinned" to the 
fluctuon energy and almost independent on n due to exponential dependence 
of the DOS (138) on E. 

5 Conclusions 

Let us resume on the main results obtained. Due to complexity of the problem 
of the electron states near quantum critical point (QCP) it is hardly believable 
that this problem can be treated rigorously. To obtain first insight into this 
we used variational approach within Feynman path integral formalism. Orig- 
inally, this approach was developed in the connection with polaron in ionic 
crystals and proved to give excellent results [24,25]. For the case of classical 
critical point (CCP) we have checked the reliability of this approach by fairly 
independent Green function method. 

The results on the electron ground state at QCP turn out to be crucially 
dependent on the anomalous space dimensionality d = D — 2 + r\ and dy- 




(152) 
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namical critical exponent z. The most interesting result is nonexistence of 
regular perturbation theory for the ground state energy for arbitrary small 
coupling constant g. In such cases singular perturbation theory emerges with 
the expansion in non-integer powers of g. For z > 2, those cases fall into range 



< d< 1. 

In the above mentioned singular perturbation-theory cases, as well as in gen- 
eral situation at large enough g (strong coupling regime) the leading term in 
the ground state energy is independent of z and is given by Eq.(97). This 
result is valid for g 2 <C Wuj (W is the electron bandwidth and u is a typ- 
ical fluctuation energy) which in fact is a criterion of consistence of contin- 
uum approximation. Physically this means that the size of self-trapped state 
(fluctuon) is much larger than interatomic distance. Otherwise a small-radius 
fluctuon likely forms, which should be considered by different methods. 

In contrast with the quantum case, at CCP the fluctuon states form a con- 
tinuum in the DOS tail. In this case the variational ffuctuon's free energy by 
Feynman method simply gives a position of the electron chemical potential in 
the tail counted from the bare band edge. The tail capacity proves ( ) 4 d 



times a numerical constant; if the electron concentration is much larger than 
this estimate the ffuctuons can scarcely contribute to the electron properties 
of material near CCP. 

References 



d + z — 2 < 1. For z < 2 it occurs at z > |(d + 2) which is consistent if 




[1] J. A. Hertz, Phys. Rev. B 14 (1976) 1165. 



42 



[2] S. Sachdev, Quantum Phase Transitions, Cambridge University Press, 
Cambridge, 1999. 

[3] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, Rev. Mod. Phys. 69 
(1997) 315. 

[4] R. B. Laughlin, G. G. Lonzarich, P. Monthoux, and D. Pines, Adv. Phys. 50 
(2001) 361. 

[5] M. Vojta, Rep. Prog. Phys. 66 (2003) 2069. 

[6] D. E. Grupp and A. M. Goldman, Science 276 (1997) 392. 

[7] A. Abanov, A. V. Chubukov, and J. Schmalian, Adv. Phys. 52 (2003) 119. 

[8] S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, J. M. Tranquada, A. 
Kapitulnik, and C. Howald, Rev. Mod. Phys. 75 (2003) 1201. 

[9] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. Bloch, Nature 415 
(2003) 39. 

[10] S. A. Grigera, R. S. Perry, A. J. Schofield, M. Chiao, S. R. Julian, G. G. 
Lonzarich, S. I. Ikeda, Y. Maeno, A. J. Millis, and A. P. Mackenzie, Science 
294 (2001) 329. 

[11] P. Gegenwart, J. Custers, C. Geibel, K. Neumaier, T. Tayama, K. Tenya, O. 
Trovarelli, and F. Steglich, Phys. Rev. Lett. 89 (2002) 056402. 

[12] J. R. Stewart, B. D. Rainford, R. S. Eccleston, and R. Cywinski, Phys. Rev. 
Lett. 89 (2002) 186403. 

[13] N. Doiron-Leyraud, I. R. Walker, L. Taillefer, M. J. Steiner, S. R. Julian, and 
G. G. Lonzarich, Nature 425 (2003) 595. 

[14] P. Coleman, C. Pepin, Q. Si, and R. Ramazashvili, J. Phys.: Condens. Matter 
13 (2001) R723. 



43 



[15] P. Coleman and A. J. Schonfeld, Nature 433 (2005) 229. 

[16] R. Blinc and B. Zeks, Soft Modes in Ferroelectrics and Antiferroelectrics, North- 
Holland, Amsterdam, 1974. 

[17] S. V. Vonsovsky and M. I. Katsnelson, Quantum Solid State Physics, Springer, 
Berlin, 1989. 

[18] G. D. Mahan, Many-Particle Physics, Plenum Press, New York, 2000. 

[19] E. L. Nagaev, Physics of Magnetic Semiconductors, Mir, Moscow, 1983; E. L. 
Nagaev, Phys. Rep. 346 (2001) 388. 

[20] M. A. Krivoglaz, Uspekhi Fiz. Nauk 111 (1973) 617 [Engl. Transl.: Sov. Phys. 
Uspekhi 16 (1973) 856]. 

[21] M. I. Auslender and M. I. Katsnelson, Teor. Matem. Fizika 43 (1980) 261 [Engl. 
Transl: Theor. Math. Phys. 43 (1980) 450]. 

[22] M. I. Auslender and M. I. Katsnelson, J. Magn. Magn. Mater. 241 (1981) 117. 

[23] M. I. Auslender and M. I. Katsnelson, Dokl. AN SSSR 265 (1982) 601 [Engl. 
Transl.: Sov. Phys.: Doklady 25 (1982) 601]. 

[24] R. P. Feynman, Phys. Rev. 97 (1955) 660. 

[25] R. P. Feynman and A. R. Hibbs, Quantum Mechanics and Path Integrals, 
McGraw Hill, New York, 1965. 

[26] M. H. Degani, O. Hipolito, R. Lobo, and G. A. Farias, J. Phys. C 19 (1986) 
2919. 

[27] W. B. da Costa and F. M. Peeters, J. Phys.: Condens. Matter 8 (1996) 2173. 

[28] L. S. Schulman, Techniques and Applications of Path Integration, Wiley, New 
York, 1981. 

44 



[29] H. Kleinert, Path Integrals in Quantum Mechanics, Statistics and Polymer 
Physics, World Scientific, Singapore, 1995. 

[30] I. S. Gradstein and I. M. Ryzhik, Table of Integrals, Series and Products, 6th 
ed., Academic Press, San Diego, 2000. 

[31] P. C. Hohenberg and B. I. Halperin, Rev. Mod. Phys. 49 (1977) 435. 

[32] A. B. Migdal, Theory of Finite Fermi Systems and Applications to Atomic 
Nuclei, John Wiley k Sons, New York, 1967. 

[33] S. Doniach and E. H. Sondheimer, Green's Functions for Solid State Physicists, 
Imperial College Press, London, 1998. 

[34] B. I. Shklovskii and A. L. Efros, Electronic Properties of Doped Semiconductors, 
Springer, Berlin, 1984. 

[35] M. I. Auslender, N. G. Bebenin, and V. P. Kalshnikov, phys. stat. sol. (b) 114 
(1982) K147. 

[36] M. I. Katsnelson and A. V. Trefilov, Z. Phys. B 80 (1990) 63. 



45 



